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Abstract. Single- wall carbon nanotubes are almost ideal systems for the investi- 
gation of exotic many-body effects due to non-Fermi liquid behavior of interacting 
electrons in one dimension. Recent theoretical and experimental results are reviewed 
with a focus on electron correlations. Starting from a microscopic lattice model we 
derive an effective phase Hamiltonian for conducting single-wall nanotubes with 
arbitrary chirality. The parameters of the Hamiltonian show very weak dependence 
on the chiral angle, which makes the low-energy physics of conducting nanotubes 
universal. The temperature-dependent resistivity and frequency-dependent optical 
conductivity of nanotubes with impurities are evaluated within the Luttinger-like 
model. Localization effects are studied. In particular, we found that intra-valley and 
inter- valley electron scattering can not coexist at low energies. Low-energy proper- 
ties of clean nanotubes are studied beyond the Luttinger liquid approximation. The 
strongest Mott-like electron instability occurs at half filling. In the Mott insulat- 
ing phase electrons at different atomic sublattices form characteristic bound states. 
The energy gaps of 0.01 — 0.1 eV occur in all modes of elementary excitations. We 
finally discuss observability of the Mott insulating phase in transport experiments. 
The accent is made on the charge transfer from external electrodes which results 
in a deviation of the electron density from half-filling. 



1 Introduction 

Single-wall carbon nanotubes (SWNTs) are cylindrical fuUerene structures 
with diameters in the nanometer range and lengths of few micrometers 
Experimental demonstration of electron transport through SWNTs ||^,^ has 
been followed by observations of atomic structure PJ^,^ , one-dimensional van 
Hove singularities [^jsj , standing electron waves and electron correlations 
[pp|,p^ in these systems. Moreover, first prototypes of SWNT functional de- 
vices - diodes and field effect transistors - have been demonstrated 
recently. 

On one hand, SWNTs can be viewed as giant macromolecules whose prop- 
erties can be learned from the first principle calculations. On the other hand, 
SWNTs are perfect one-dimensional (ID) model systems to be studied by 
methods of the solid state theory. This somewhat reductionist but insightful 
approach provides a reasonable description for the bulk of experimental data 
obtained up to date. 
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On a single-particle level physical properties of SWNTs are determined 
by their geometry. Depending on the wrapping vector, SWNTs can either be 
ID metals or semiconductors with the energy gap in sub-electronvolt range. 
This has been confirmed by direct observation of ID van Hove singularities 
in scanning tunneling microscopy experiments [^j^] . 

In this paper we address the role of the Coulomb interaction in ID SWNTs 
and review recent results on electron correlation effects. Away from half-filling 
electron correlations are well described by Luttinger-like models ^3 
particular, the non-Fermi liquid ground state of the system is characterized 
by a power-law suppression of the density of electronic states near the Fermi 
level. This effect has been observed in single- and, presumably, multi-wall 
[H nanotubes, as well as injunctions between metallic SWNTs Transport 
properties of metallic nanotubes with impurities are affected by the Coulomb 
interaction as well |p^ . 

The low-energy properties of SWNTs are different at half-filling due to 
the umklapp scattering. The latter is coupled to the strongly interacting 
total charge mode. This makes the umklapp scattering a strongly relevant 
perturbation. As a result, the Mott-like electron instability occurs in the 
electronic spectrum of SWNTs and energy gaps open in all modes of the 
elementary excitations [^^. The transition to the Mott insulating phase 
should manifest itself by an increased resistivity of half-filled SWNTs at low 
temperatures. 

The plan of the paper is as follows. In Section 2 we derive an effective low- 
energy model for metallic SWNTs with arbitrary chirality, which naturally 
includes one-electron backscattering on impurities as well as two-electron 
umklapp scattering at half-filling. We discuss the non-Fermi liquid properties 
of clean SWNTs within the Luttinger-like model in Section 3 and evalu- 
ate electron transport in SWNTs with impurities in Section 4. The effect of 
interactions beyond the Luttinger model is analyzed within the renormaliza- 
tion group scheme in Section 5. The Mott-like electron instability is further 
investigated in Section 6 using the self-consistent harmonic approximation. 
Finally, we discuss observability of the Mott-insulating phase in experiments 
(Section 7). 
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2 Universal Model of Metallic Nanotubes 



2.1 Microscopic theory 

Structurally uniform SWNTs can be characterized by the wrapping vector 
w = N^a^ + N^a^ given by the linear combination of primitive lattice 
vectors a± = (±1, \/3)a/2, with a ~ 0.246 nm (Fig. |^). It is natural to sep- 
arate non-chiral armchair {N^ — N-) and zig-zag {N+ — —N^) nanotubes 
from their chiral counterparts. Recent scanning tunneling microscopy studies 
[^H^ have revealed that individual SWNTs are generally chiral. According 
to the single-particle model, the nanotubes with iV+ — N- = mod 3 have 



Fig. 1. Graphite lattice consists of two atomic sublattices p = denoted by 

filled and open circles. SWNT at the angle x to a; axis can be formed by wrapping 
the graphite sheet along w = OO' vector. 

gapless energy spectrum and are therefore conducting; otherwise, the energy 
spectrum is gapped and SWNTs are insulating. 

We consider metallic SWNT whose axis x' forms an angle x = arctan[(A^_ — 
N+)/^/3{N+ +iV_)] with the direction of the chains of carbon atoms {x axis 
in Fig. |l|). We expand the standard single-particle Hamiltonian ||l8| Hk for 
electrons on two atomic sublattices p = ± of a graphite sheet (Fig. |l|) near 
the Fermi points aK (with the valley index a = ±, and K = (47r/3a,0)) to 
the lowest order in q = A; — aK = (7(cosx,sinx). Introducing slowly varying 
Fermi fields ippas{x') = L^^/^ J2q=27rn/L e^'^'^' apas{q + aK), we obtain, 

Hk = -iw^ae^'f"^ Ux'i^l^,d,,i>-po.s, (1) 

pas 

where v ~ 8.1 x 10^ m/s is the Fermi velocity, and s = ± is the electron spin. 
The kinetic term can be diagonalized by the unitary transformation 

to the basis (pras of right- (r = and left- {r = —) and moving electrons. 
The Coulomb interaction has the form 

Htnt^^ Vpp-{2dK) J dx'^pl^^,^l,^^^,^p'a-,s'^pais, (3) 

pp' ^{ai} ,ss' 

with the matrix elements Vpp/ {2aK) corresponding to the amplitudes of intra- 
{p = p') ctnd inter- (p = —p') sublattice intra- (a = 0) and inter- (a = ±1) 
valley scattering (here a — (ai — ai)/2 = {a^ — a2)/2). In Eq. (||) we assume 
that the fields ip are varying slowly on the scale of the screening radius Rs of 
the Coulomb interaction determined by the distance to a gate electrode. This 
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Table 1. Scattering amplitudes AV{0), Vpp{2K), Vp-p{2K) in units ae^ /2iikR for 
all SWNTs with 2R/a = 4 — 7. The cutoff parameter ao = 0.526a is obtained from 
the requirement that the on-site interaction in the original tight-binding model is 
equal to the difference between the ionization potential and electron affinity of sp^ 
hybridized carbon |l9[. 



ao/a 


AV{0) 


Vpp{2K) 


\Vp-p{2K)\ 


0.4 
0.526 
0.7 


0.44265 - 0.44274 
0.17378 - 0.17395 
0.04880 - 0.04895 


0.97060 - 0.97095 
0.53549 - 0.53561 
0.24778 - 0.24797 


0.6 - 2.2 X 10"^ 
0.5 - 1.6 X 10"^ 
0.3-1.5 X 10"^ 



corresponds to a contact-interaction approximation. Equation (^) is therefore 
valid at large length scale x' ^ Rg and low electron energy E ^ Tiv/Rs- 

The dominant contribution to the intra- valley scattering amplitudes Vpp' (0) 
comes from the long range component of the Coulomb interaction, Vpp(0) ~ 
Vp-p{0) ~ eVC ~ {2e^/K) ln(i?^/i?), C being the capacitance of SWNT per 
unit length The differential part AV{0) = Vpp{0) - Vp-p{0) of intra- 
valley scattering as well as the intra-sublattice inter- valley Vpp{2K) are esti- 
mated at AV{0),Vpp{2K) - ae^/nR . Despite AV [Q) ,Vpp{2K) being much 
smaller than Vpp{Q), they cause non-Luttinger terms in the low-energy Hamil- 
tonian which will be important in the further analysis. The matrix elements 
Vpp'{2aK) have been evaluated numerically for all SWNTs with radii R in 
the range 2R/a = 4-7 {2R/a = 5.5 for (10,10) SWNTs). We found that 
dimensionless amplitudes 2ttkR[AV{0), Vpp{2K)]/ae^ show very weak depen- 
dence on the radius of SWNT and its chiral angle (see Table 1). The results 
are sensitive to the value the short distance cutoff ao ~ a of the Coulomb 
interaction. 

The inter-sublattice inter- valley scattering amplitude Vp-p{2K) is almost 
three orders of magnitude smaller than AV{Q), Vpp{2K). This is due to the 
C3 symmetry of a graphite lattice {Vp^p{2K) = for a plane graphite sheet). 
The matrix elements Vp-p{2K) are generally complex due to the asymmetry 
of effective ID inter-sublattice interaction potential (the matrix elements are 
real for symmetric zig-zag and armchair SWNTs). Let us note that after the 
unitary transformation (||) of the Hamiltonian H = Hk + i?mt, the chiral 
angle x enters only into the inter-sublattice inter-valley scattering matrix 
elements Vp-p{2K) . Due to the smallness of these matrix elements, the low- 
energy properties of chiral SWNTs are expected to be virtually independent 
of the chiral angle. 
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2.2 Bosonization 



The transformed Hamiltonian H can be bosonized by introducing the phase 
representation of the Fermi fields , 



irqpx' + — {Oas + r4ias} 



(4) 



The phase variables 9asi 4'as are further decomposed into symmetric 5 — + 
and antisymmetric 5 = — modes of the charge p and spin a excitations, 
Oas = Op+ + sOcr+ + ctOp- + asOa- ,0 — 9,(1). The bosonic fields satisfy the 
commutation relation, [0jsixi),(l)j>s'{x2)] — i(7r/2)sign(a;i — X2)Sjj'Sss' ■ The 
Majorana fermions rjras are introduced to ensure correct anticommutation 
rules for different species r,a,s of electrons, and satisfy ['r]ras,'nr'a's']+ = 
6ss'- The quantity qp = 7rn/4 is related to the deviation n of the 
average electron density from half-filling, and a ~ a is the parameter of the 
exponential ultraviolet cutoff. 

Neglecting the inter-sublattice inter- valley scattering we arrive at the uni- 
versal phase Hamiltonian of metallic SWNTs, 

^ = E Elf / dx'{K-l{d.,e,sf+K,s{d.'C^,s?} + ^ f dx' 

{[AV{Q) ~ Vpp{2K)][cos{AqFx' + 29 p+) cos29„+ - cos29p^ cos26'^_] 
-AV{0) cosiAqpx' + 29 p+) cos 29 p^ + AV{<d) cos{AqFx' + 29 p+) cos 29^^ 
-AV{Q) cos 29„+ cos 29 p- + AV{Q) cos 29„+ cos 26*^- 
-Vpp{2K) cos{AqFx' + 29p+) cos 2(f>a- + Vpp{2K) cos 29„+ cos 2(/)„_ 
-HFpp(2X)cos26'p_ cos20^_ +Vpp{2K)cos29„-Cos2(i)„^}, (5) 



Vjs = AjsBjs and Kjg — ^ Bjs /Ajg being the velocities and interaction 
parameters for different modes j, 6 of excitations. The parameters Ajs, Bjs are 
given by Ap+ = 1 + AV{0)/2~Vpp{2K)]/2nv, A^s = 1 - [Z\F(0)/2-h 

5Vpp{2K)]/2T:v, B,s = 1 + AV{Q)/4nv, with V{0) = [Vpp{0) + Vp^p{0)]/2. 
The renormalization of the parameters Kjs, v^s by the Coulomb interaction 
is strongest in p+ mode. Assuming k — 1.4 R ~ 0.7 nm, and Rs — 100 
nm we obtain Kp^ ~ 0.2. The interaction in the other modes is weak: Kjs = 
l + 0{a/R). 

Let us note that the Hamiltonian (^) has the same form as the phase 
Hamiltonian for a two- leg Hubbard- type ladder provided that the dif- 
ference in definitions of the fields and (j = p, a) in terms of densities 
of right- and left-movers in two energy bands is taken into account. 



2.3 Impurity scattering 

Disorder in the atomic potential is described by the Hamiltonian pl| , p^ , 
Himp = J2paa'sl dxVpa{x)'4>lasi'pa' s ■ Here the impurity potential Vp{x) at 
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the sublattice p = ± is decomposed into intra-valley {a = {a' — a)/2 = 0) 
and inter-valley (a = ±1) scattering components. First, we transform the 
Hamiltonian Himp to the basis of right- and left-movers (r = ±), see Eq. 
(|). When the range of the impurity potential is much larger than the lattice 
constant, the backward scattering (r — r) is ineffective [Q. We consider 
the case of short range impurity potential and retain backscattering terms 
in the Hamiltonian. The forward scattering is discarded because it does not 
contribute to the transport. 

In the limit of weak impurity potential, the interaction between the elec- 
trons and the impurities can be parameterized by uncorrelated Gaussian ran- 
dom fields, 'r]{x) and ^(x) expressing the intra-valley and the inter- valley back- 
ward scattering, respectively. The fields satisfy {r]{x)ri{x'))imp = Di6{x — x') 
and {i,(x)£,* {x')) imp = D25{x — x'), where (• • ■)imp is the configurational av- 
erage. The factors Di and D2 are given by u/ti and v/t2 with the scattering 
time Ti (T2) due to the intra-valley (inter- valley) backscattering. The Hamil- 
tonian of impurities is given by 



Hirap 



j dx-q{x) ^ ipl^^lp^ras + j dxl ^{x) ^ Vj+>-r-s + h.C. \ . (6) 
ras \ rs ) 



Note that the intra-valley (inter- valley) backward scattering is parameterized 

limp = -^imp + Hfmp 



by a real r]{x) (complex ^{x)) field. The Hamiltonian Hi,np = H} + 



is expressed in terms of the phase variables as follows, 
Hi,., = ||/<i.„.)^™exp,-2i.„., 

ras 

X exp {-\r{6p+ + s6„+ + a9p_ + as9a-)} , (7) 
" / da; ^ r exp {-ir{2qFX + Op^ + 89^+)} 

rs 

X [^{x)exp{-i{(f>p_+s(l)„_)} + h.c.]. (8) 

with the Pauli matrices ctj,, <7z originating from the Majorana representation 
of the Fermi fields (0). 



3 Luttinger model limit 

The phase Hamiltonian (||) of clean metallic SWNT consists of a part quadratic 
in bosonic fields describing scattering of electrons within the same branch of 
the spectrum and the non-quadratic part describing inter-branch scatter- 
ing. In the Luttinger model limit one neglects the differential part AV{0) of 
intra-valley scattering and the inter- valley scattering Vpp'{2K). In this case 
the Hamiltonian (^) reads, 

= E E If / {Kis'id^'S.sr + K,sid.,q^,sr} , (9) 
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with parameters vjs = Vy/Ajs and Kjs = 1/ y/Ajs determined by Ap+ = 
1 + 4y{0)/TTv, Ajs = 1 for {jS) ^ {p+). The parts of the Hamiltonian (|) 
in the four sectors of excitations {j6) are decoupled. In particular, using the 
relations p = (2/7r)9x^p+, j — {2v/Tr)dx4>p+ for the charge density ep and the 
electric current ej (e > 0) in SWNT, the Hamiltonian (||) in the total charge 
p+ sector can be rewritten as follows, 



Hp+ = - I dx' 



1 2 1 n 

2^ +2 



dx'dx"p{x')Vix' - x")pix"), (10) 



(in Eq. (^ the contact interaction approximation V{x) = V{0)S{x) has been 
made). In Eq. (|l^) j/ = i/nv is the density of electronic states in SWNT; the 
first term describes the kinetic energy of right and left moving electrons and 
the second corresponds to the Hartree part of the interaction. 

Electronic properties of SWNTs in the Luttinger limit have been inves- 
tigated in e.g. Refs. [|l3yi4|] . The density of electronic states near the Fermi 
level is suppressed in a power-law fashion. This is a signature of the orthog- 
onality catastrophe which occurs due to the non-fermionic (bosonic) nature 
of low-energy excitations in SWNTs. As a result, the differential conduc- 
tance in the tunneling regime shows a power-law dependence on tempera- 
ture, dl/dV (X T", for < T and voltage, dl/dV cx for eV > T. The 
exponent a is given by a = i-^p+ + ^'^p+ ^ 2)/8 for the tunneling from a 
metallic electrode into the "bulk" of SWNT, a = {Kp_l - l)/4 for the tunnel- 
ing into the end of SWNT, and a — {Kp^ — l)/2 for the end-to-end tunnehng 



in nanotube heterojunctions |10 . Moreover, scaled differential conductance, 
{dI/dV)/T°', is a universal function of the ratio eV/T . The collapse of the 
data taken at different temperatures to a single universal curve provides a 
comprehensive check of validity of the Luttinger model for SWNTs j|Jl^ . 



4 Effect of impurities 

In this section, the transport properties of metallic SWNTs with impurities 
are studied in the Luttinger model limit p6| . The dynamical conductivity 
17(0;) is expressed via the memory function AI{lj) as follows |2^ , 

^ ' uj + Miu;)' ^ ' -x(0) ^ ^ 

where = -ijdx dte(''^~'')*([A(a;, i), A(0, 0)]) with 77 -> +0, (• • •) 

denotes the thermal average with respect to the Hamiltonian H ~ HL+Himp, 
F = [j,H] with j being the current operator, and x(0) = {j',j)ui=o — — 4w/7r. 
The temperature dependence of the resistivity is given by, 

rHiKp++3)/4) /27rTy^^+-^)/^ 

^ = ^-r((x,, + 3)/2) [—) ' (12) 
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where pso = 2(^'''*)^^' ^^'^ resistivity of the non-interacting 

system in the Born approximation, lof = v/d, and r(z) is the gamma func- 
tion. It is remarkable that p/pso is independent of the scattering strength. 
The repulsive Coulomb interaction (A'p+ < 1) leads to an enhancement of 
the resistivity at low temperatures. For typical nanotubes with Kp^ ~ 0.2, 
the resistivity scales as p cx T~°-^. 

The optical conductivity Kealuj) behaves like at high frequen- 

cies w ^ T and as at low frequencies. The position of the peak 

in Recr(w) is given by. 



uJFTi T2' r{{Kp+ + 3)/2) 



(13) 



We will further study localization of electrons which corresponds to the 
pinning of the phase fields 9js,4'jS- The localization is described within the 
renormalization group (RG) formalism. RG equations can be derived by as- 
suming scale invariance of correlation functions along the lines of Ref . [ p3| , 

(ViY = {3 - iKp+ + + Kp^ + K,^)/2}Vi, (14) 

W - {3 - {Kp+ + + K;^ + K-^)/2] V2, (15) 

{Kj+y = - (Pi/Xi + V2/X2) K^+Uj+, (16) 

{uj+y = - (Pi/Xi + V2/X2) K,+u]+, (17) 

(if,_)' = - {ViK]_/Xi - V2/X2) uj^, (18) 

(u,_)' = - {ViKj^/Xi + V2KJ}/X2) u]_, (19) 

where ( )' denotes d/d£ with M = dln(a/o) (a is the new lattice constant), 

^1 — Wp+ U^j^ Up_ , VL2 — Up+ "p- "(T- 

j = p OT a. The initial conditions for the above RG equations are as fol- 
lows, V,{0) = D,d/{7rv^), Kp+{0) = u'l = Kp+, and K,+ {0) - Kp^{0) - 
K,-{0) = u,+ {0) = up_(0) = M,_(0) = 1. 

The temperature dependence of the resistivity. Fig. ^, can be evaluated 
by solving the RG equations. The enhancement of the resistivity at low- 
temperatures is due to the electron localization. Solution of RG equations 
shows that the intra- valley (inter- valley) scattering pins the phases, 0p+, 9a-+-, 
9p-, and 0^^ {Op+, Oa-+, (t^p-, and ipu-)- Since the conjugate variables, 0p- 
and (/>p-, or and (j)^- cannot be pinned at the same time, the localization 
due to two kinds of the scattering cannot occur simultaneously. For the pa- 
rameters of Fig. ^, the intra-valley scattering dominates over the inter-valley 
scattering at low temperatures (see inset of Fig. ^). 

The high frequency behavior of the optical conductivity is not modified 
by the effects of the localization. The power law behavior at low frequencies 
survives due to a finite density of states at the Fermi energy (similarly to the 
non-interacting case pJ). 




Fig. 2. Temperature dependence of the resistivity p of SWNT witli the following 
parameters, hujfti = 300 and ■nu}FT2 ~ 600. The solid (dotted) line corresponds 
to RG analysis (perturbation theory). Inset : The ratio pij p\ as a function of 
temperature. 



5 Effect of interactions beyond the Luttinger model 

Until now the Coulomb interaction in SWNTs has been treated within the 
Luttinger model. In order to describe the interaction effects beyond this ap- 
proximation, one has to consider the full Hamiltonian (|^). The low energy 
properties of Eq. (^ can be investigated by the RG method (see Rcf. for 
details). At half-filling, qp = 0, we obtain the following RG equations. 









(20) 


(^.+ )' 


= -{Kl^myl+yl+yl- 


^y|) , 


(21) 




= -{Kl^l8){yl + yl + yl- 


^yg) , 


(22) 




= -{Kl_myl + yl + yl) 


+ (l/8)(2/2+y|+y2) ^ 


(23) 


(yiY 


= {2- {Kp+ + 


K„+)}yi - 


{y2y6 + yayb + yvys) / 4: , 


(24) 


{V2)' 


= {2- {Kp+ + 


Kp-)}y2 - 


(yiye + y3y4 + yryd)/^ , 


(25) 




= {2- {Kp+ + 


Ka-)}yz - 


(yiys + y2y4)/4 , 


(26) 


(y^y 


= {2- {Kp^ + 


K.-)}yi- 


(y2y3 + y5y6)/4 , 


(27) 


(ys)' 


= {2- {K,+ + 


Ka-)}y5- 


(yiys + y4y6)/4 , 


(28) 


(ye)' 


= {2 - {K„+ + 


Kp-)}y&- 


(yiy2 + y4y5 + y8y9)/4 , 


(29) 


(yr)' 


= {2- {Kp+ + 


l/^.-)}y7 


- (yiys + y2y9)/4 , 


(30) 


(ys)' 


= {2 - (A',+ + 


i/^.-)}y8 


- (yiy? + y6y9)/4 , 


(31) 


(yg)' 


= {2 - (i^p_ + 


l/K.-)}y^ 


- {y2y7 + yeys) / 4: . 


(32) 
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The initial conditions for Eqs. (|o|)-(||) are Kjs{0) = Kjs, yi = [^V{0) - 
Vppi2Ko)]/{7Tv), V2 = -ys - -2/5 = 2/6 = - AV {0) / {irv) , y^ = [Vpj,i2K) - 
AV{0)]/{ttv), yr = -2/8 = ~Vpp{2K) / {nv) , and yg = Vppi2K) / (ttv) . In de- 
riving the RG equations, the non-hnear term cos 20(j_ cos 20o-_ is omitted 
because this operator stays exactly marginal in all orders and is thus decou- 
pled from the problem p^. The RG equations away from half-filling can be 
obtained from Eqs. (pO[)-(|32D by putting yi, y2, 2/3 and 2/7 to zero. Hereafter 
we concentrate on the case N = IQ, k = 1.4, Rg = 100 nm and ao/a = 0.526 
and estimate the initial values of the RG parameters using Table 1. 

Away from half-filling, the quantities ifo-+, -R'p-, and K~}i renormal- 
ize to zero and the coefficient of cos 29^+ cos 26 (cos 26*0-+ cos 20o-- and 
cos 20p_ cos 200— ) tends to —00 (00). As a result, the phases do-+,6p- and 
are pinned at (0^^+ , , 0o-- ) = (0, 0,7r/2) or (7r/2, 7r/2, 0) so that the 
modes cnt and p— are gapped. In this case, the asymptotic behavior of the 
correlation functions at x ^ 00 is determined by the correlations of the gap- 
less p+ mode, (^c'"'^P+(^)e-''^Op+(°'>) - 2;-'^'^p+/2 and (^e'"^'i>P+i^)e-'"'^p+(°')) 

^-m jiKp^ (n~\ and 2 correspond to 2qp and ^qp density waves and m = 1 
to a superconducting state). Since -ft'p+ « 0.2, the 2qp density wave correla- 
tions seem to be dominant. However, we found that the correlation functions 
of any 2qp density wave decay exponentially at large distances due to the 
gapped modes. We therefore are looking for the four-particle correlations. 
The 4gi? density waves dominate over superconductivity for ifp+ < 1/2 
[ p5[|2^ . Such density wave states are given by the product of the charge 
p^(x)p-(x) or spin Sj^[x)S-(x) densities at different sublattices. Substitut- 
ing the values of the pinned phases we observe that pj^{x)p^(x) vanishes, 
and the dominant state is the \qp spin density wave with correlation func- 
tion (5+(a;)S'_(a;)5+(G)S'_(0)) - co^ \qpx j x'^^ •'^ . 

At half-filling the solution of the RG equations (^)-(^) indicates that 
the phase variables ^p+, ^0-+, 6p-^ and are pinned and the all kinds of 
excitation are gapped. In other words, the ground state of the half-filled AN 
is a Mott insulator with spin gap. The same conclusion has been drawn from 
the model with short range interactions ||2^,^ . The pinned phases are given 
by (6'p+,6'<,+ ,6lp_, (0,0,0,0) or (7r/2, 7r/2, 7r/2, 7r/2) since the first, 

second and 6-9-th non- linear terms in Eq. (^) scale to —00. The averages 
{p+{x)p-{x)) and {S+{x)S-(x)) are both finite, which indicates formation 
of bound states of electrons at different atomic sublattices. 

The states derived from the present analysis are characteristic for the 
long range Coulomb interaction. In fact, for the on-site plus nearest neighbor 
interaction the dominant states correspond to the density waves at half-filling 
and to the superconducting state or the density waves away from half-filling 
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6 Mott-insulating phase 

To estimate the gaps in different modes of excitations at half-filling, we em- 
ploy the self-consistent harmonic approximation which follows from Feyn- 
man's variational principle We consider a trial harmonic Hamiltonian of 
the form: 

^0 = E /d^' { Kjsiid.'d.sY + (1 - 5,Js-)q%e%] 



+ K^sm,cj)^sf + S,Js-ql^l]}, (33) 



QjS being variational parameters. By minimizing the upper estimate for the 
free energy F* = Fq + {H — Hq)q one obtains the following self-consistent 
equations, 

ql+ - ^-f^^ Cp+[Vpp{2K) - Z\F(0)]c,+ + Z\F(0)cp_ + Vpp{2K)da- {M) 

IT Oj U p-)- 

2K 

ll- = —rj^Cp^^V{0)cp+ + AV{Q)c„+ - Vpp{2K)d^_, (35) 
ql+ = c.+ [Vpp{2K) ~ AVmcp+ + Z\F(0)cp_ - Vpp{2K)d,-{?,Q) 

IT Qi V fj 

2 

ql- = —rjTz da-Vpp{2K) {cp+ - - Cp_} , (37) 

wherecji = {cos29js)„ = cos20'^'g\-fdqjs)'^'\ d^y- = (cos 20cr-)n_= cos2(f>''J'l' (jaqa-)'^^'^' 
(...)q denotes averaging with respect to the trial Hamiltonian (|33|), and 7 ~ 
0.890. Note that (cos20o.-)q — 0, so that only the terms of the Hamiltonian 
(||) which scale to the strong coupling (see end of Section ^) contribute to 
Eqs. (H-dl). 

In the limiting case of interest, |ZiF(0)|, |V^p(2i^)| ^ v and Kp^ <C 1, the 
solution of Eqs. (^4|)-(^7|) can be found in a closed form, giving rise to the 
following estimates for the gaps Ajs = Vjsqjs in the energy spectra. 







(38) 




7a V Trvp+ J 


^p- = 


Vp+ 


(39) 


A^+ = 


\Vpp{2Ko)-AVm 
Vp+ 


(40) 


A,- = 


\Vpp{2K,)\ ^ 


(41) 



withl/p+ = {[AV{0) - Vpp{2K)f + [AV{Q)Y +[Vpp{2K)YY'^ . In the above 
expressions we used the approximation, v/vp+ — Kp^ and v/vjs = Kjs = 1 
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Fig. 3. The energy gaps Ajs for the modes p+, a—, a+, p— at ao = 0.526a (lines 
marked by crosses, from top to bottom) and for the mode p+ at oq = 0.4a (triangles) 
and at ao = 0.7a (squares). The energy is in units hv/a ~ 2.16 eV for d — a. 

for cr± and p— modes. The formulae (|3^)-(|4l|) indicate that the largest gap 
occurs in the p+ mode, albeit all four gaps are of the same order for realistic 
values of the matrix elements (see Table 1). The gaps decrease as Ajs oc 
(l/i?)^/(i-^''+) ~ with the nanotube radius. 

In Figure ^ we present numerical results for the gaps Ajg for the short 
distance cutoff of the Coulomb interaction ao = 0.526a. The data for some- 
what larger and somewhat smaller values of ao indicate possible variations of 
the gap Apjf. due to the uncertainty in the cutoff. The gaps can be roughly 
estimated at Ajs ^ 0.01 - 0.1 eV for typical SWNTs with R ~ 0.7 nm. 

The resistivity p of metallic SWNTs increases exponentially p oc exp(Z\p+ /T) 
at low temperatures T <^ Ap+. On the other hand, at high temperatures, 
T ^ Ap+, perturbation theory with respect to the non- linear terms of the 
Hamiltonian (||) gives a power-law behavior of the resistivity p ~ T'^^''+~^ /N"^ 
due to umklapp scattering at half-filling. Note that the power-law is different 
from that governing the impurity scattering, Eq. (p^. A resonant increase 
of the resistivity of SWNTs at half-filling is a characteristic signature of the 
Mott insulating phase. 

Due to the gaps in the spectrum of bosonic elementary excitations, the 
electronic density of states should disappear in the subgap region and display 
features at the gap frequencies and their harmonics. These signatures should 
be observable by means of the tunneling spectroscopy. 

7 Observability of the Mott-insulating phase 

In this section we discuss various factors which might influence the observabil- 
ity of the Mott-insulating phase in realistic systems. An exponential increase 
of the resistivity of half-filled metallic SWNTs at low temperatures can also 
occur due to deformation-induced gaps Ad in the single-particle spectrum 
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[|0|. These gaps are in the range Ad < 0.02 cV for typical SWNTs with 
R ~ 0.7 nm. This value is comparable with (but somewhat smaller than) 
our estimate for the Mott gaps. We therefore cannot exclude the interplay 
between the two mechanisms suppressing electron transport at half- filling. 
Note that the deformation induced gaps decrease as 1/i?^ with increasing 
nanotube radius which should be contrasted to dependence for 

the Mott gaps. Moreover, deformation-induced gaps strongly depend on the 
chiral angle of SWNTs, whereas Mott gaps are determined by the radius of 
SWNTs and depend weakly on their chirality. 

Isolated metallic SWNTs are half-filled systems. However, in generic ex- 
perimental layouts the nanotubes are contacted to metallic electrodes. The 
difference AW in the work functions of a metal (typically Au or Pt) and 
nanotube results in a charge transfer between the nanotube and electrodes, 
which shifts the Fermi level of the nanotube downwards away from half-filling 
(for AW > 0). 

To be specific, consider a metallic SWNT surrounded by a coaxial cylin- 
drical gate electrode of radius Kg ^ R- The nanotube {x > 0) contacts the 
2/z-plane of metallic electrode (x < 0) at x = (see inset of Fig. The in- 
teraction term of the Hamiltonian (^0|) is modified by taking the electrostatic 
potential Vg of the gate as well as the image charge at the metallic electrode 
into account. 



with the modified Bessel functions Iq, Kq. The kernel Vq describes the long- 
range Coulomb interaction, V{x) — 1/K|a;|, for i? <C |a;| <C Rs- The interac- 
tion is screened at large distances ^ Rs, so that Vq=o — i^e^ln) \n{Rs/R), 
in agreement with our previous estimate (see text below Eq. (y)). 

Minimization of the total energy given by the Hamiltonian (|l^ with the 
interaction term (^) determines the profile of the charge density ep{x) and 
the deviation —AE{x) of the Fermi level from half-filling, AE = p/v (see 




(42) 



The Fourier images of the kernels V{x), M{x) are given by, 




(43) 



Refs. |3^,|| for details). 



AE{x) 



AW \n{Rs/R) 
vV{Q) \tl{z/R) 



vV{Q) Rs 



ceVg X 



, for i? < X < i?, 



'S ; 



(44) 



AE{x) 




Y", for i?s <C a; 



(45) 
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Fig. 4. The energy AE of the gapless point (half filling) with respect to the Fermi 
level (units of AW / {2e'^ ly / k)) . The curves from top to bottom at x/R = 1000 
correspond to eVg/AW — —2,0,1,3. The screening radius Rs/R = 75, 300 for 
dashed and solid curves corresponds approximately to 50 nm and 200 nm for (10,10) 
SWNTs. The Coulomb interaction is characterized by vV (fi) / \n{Rs / R) = 5. Inset: 
layout of the system. 

where c — (l/vr) / dx\lQ{x)\~'^ ~ 1.33 and the Coulomb interaction is sup- 
posed to be strong, vV{Qi) ^ \n{Rs/R). 

Equation (^) shows that the density of charge transferred to SWNT due 
to the mismatch AW of the work- functions decays logarithmically slow with 
the distance x from the metallic electrode This is related to the poor 
screening of the Coulomb interaction in ID SWNTs. The influence of a gate 
voltage on the charge density in SWNT is suppressed by a factor x/Rs in 
the vicinity of electrode, x < Rs, cf. Eqs. (^), (^5|). In order to achieve half- 
filling condition AE = in samples with small distances d < Rs between the 
electrodes, one should compensate for this effect by applying a higher gate 
voltage. 

The coordinate dependence of the charge density (or the deviation AE{x) 
from half-filling) is shown in Fig. ^ for several values of the gate voltage and 
the screening radius. If a substantial positive voltage is applied to the gate, 
the charge density in SWNT can change sign (from positive to negative) with 
increasing distance from the metallic electrode. Let us note that the Mott 
gap should result in the formation of half-filled incompressible regions near 
such "charge neutrality points", AE{x) = 0. This is in contrast to isolated 
SWNTs where the Mott transition occurs uniformly in the whole system. 
The electron transport through the incompressible regions seems to be an 
interesting problem for future research. 
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8 Conclusions 

We have investigated manifestations of electron correlations in metallic single- 
wall carbon nanotubes. Effective low-energy model of metallic SWNTs with 
arbitrary chirality is developed starting from microscopic considerations. The 
unrenormalized parameters of the model show very weak dependence on 
the chiral angle, which makes the low-energy properties of metallic SWNTs 
chirality-indcpcndcnt. In the absence of inter-branch electron scattering the 
model corresponds to two-channel Luttinger liquid. The impurity scattering 
in such Luttinger liquid is investigated using memory function technique and 
the rcnormalization group method. We show that the intra-vallcy and intcr- 
valley backscattering can not coexist at low energies. The properties of clean 
SWNTs are investigated beyond the Luttinger liquid limit. The ground state 
away from half-filling is found to be 4qp spin density wave. At half-filling 
the umklapp scattering affects strongly interacting mode of total charge ex- 
citations. This makes the umklapp scattering strongly relevant perturbation 
which gives rise to the Mott-insulating ground state. The latter is charac- 
terized by " binding" of electrons at two atomic sublattices of graphite. The 
energy gaps in all modes of excitations are evaluated within self-consistent 
harmonic approximation. Finally, observability of the Mott phase in realistic 
SWNTs is discussed. 
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